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Abstract 

Parametric down-conversion (PDC) offers the possibility to control the fabrication of non-Gaussian 
states such as Fock states. However, in conventional PDC sources energy and momentum conserva- 
tion introduce strict frequency and photon number correlations, which impact the fidelity of the pre- 
pared state. In our work we optimize the preparation of single-photon Fock states from the emission 
of waveguided PDC via spectral filtering. We study the effect of correlations via photon number re- 
solving detection and quantum interference. Our measurements show how the reduction of mixedness 
due to filtering can be evaluated. Interfering the prepared photon with a coherent state we establish 
an experimentally measured fidelity of the produced target state of 78%. 



1 Introduction 

A practical single photon source should meet the demands of high fidelity, brightness and efficiency. 
As a source of photon pairs, parametric down-conversion (PDC) enables conditional state preparation — 
heralding on one of the twin photons (idler) gives exact timing information about the other (signal). In 
addition to this, single pass PDC sources have the advantage of operational simplicity and the quality of 
the quantum state preparation is mainly limited by the available detectors. 

Recent experiments have demonstrated a full tomographic characterization of Fock states or photon 
added/subtracted states prepared from the PDC emission in the low gain regime [T}|5). Nevertheless, 
since the photon-number distribution of signal and idler is thermal, higher photon-number components 
degrade the fidelity of the prepared target state. This effect becomes prominent at higher gains i.e. if high 
preparation rates are needed. In addition to photon-number correlation, spectral correlation between the 
PDC photon pairs cannot be avoided in conventional sources ||6jj9|. Although the spectral correlation is 
a useful resource for the generation of highly entangled states flQf[L2) , the preparation of non-Gaussian 
states in one of the twin modes leads to mixedness whenever the other field mode is traced out | [L3| . 

In the past, the manipulation of the spectral properties has been in the interest of broadband and 
efficient non-linear classical frequency conversion fl4). More recently, tailoring these properties has 



been suggested as a method of source engineering for quantum optics 1 15 16 1. Several sophisticated 
methods have been proposed for state decorrelation, which unfortunately often require means that are 
available only under specific conditions p7] - |2T| . Nevertheless, the decorrelation of the PDC photons 
is always possible via spectral filtering, which drastically changes the shape of the spectral correlation 
function. The influence of the spectral correlation on the purity of a single photon state prepared by 
filtering was first experimentally demonstrated by Wasilewsky et al. J22) . 

In addition to the optimal state preparation, a precise characterization of the prepared photonic states 
both in photon-number and spectral degrees of freedom is required. The latest developments in pho- 
ton counting technologies allows to characterize photonic states via the direct detection of the photon- 
number statistics. One of these techniques applicable for pulsed light is time-multiplexed detection 
(TMD), which enables a loss-tolerant reconstruction of photon statistics [23,24]. Regarding the spectral 
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properties, their impact on the state preparation can be studied via Hong-Ou-Mandel (HOM) type inter- 
ference, originally applied to observe the temporal length of PDC photons [25|. This experiment can be 
implemented with a symmetric beam splitter and a coincidence counting device. Two indistinguishable 
photons impinging simultaneously on the different ports of the beam splitter bunch together and no co- 
incidences are detected. Thus, the measurement of HOM interference dip between signal and idler is a 
measure of their spectral indistinguishability, or in other words of their equity [12j. However, for char- 
acterizing the fidelity of the prepared single photon, one must employ an independent, spectrally pure 
reference field, which can be a highly attenuated coherent state. The measurement of HOM interference 
between these two independent sources allows to quantify their spectral overlap, as evidenced in several 
studies |26]-p0|. Considering the direct probing of more evolved non-Gaussian states, such as displaced 



single photons, the determination of this overlap is essential 1 31 1. 

Lately, PDC sources with extreme high brightness have been accomplished by employing waveg- 
uided structures rather than bulk crystals [ 32 1 . As filtering happens at the cost of photon flux, waveguides 
(WGs) can offer benefits for preparing high fidelity single photons. In this paper we show how to opti- 
mally adapt the bandwidth of the filter such that the trade-off between the fidelity and brightness of the 
source is balanced. In Sec. [2] we provide easy applicable tools for analyzing and estimating the spectral 
correlation of a waveguided PDC source and give a recipe of how to design a filter for optimal state 
decorrelation. The spectral properties and the performance of the employed filters are verified in Sec. [3] 
Our main focus, the characteristics of the filtered state is presented in Sections |4] and [5j we study the 
statistics of the heralded state and finally determine the spectral overlap between signal and reference. 
This allows us, to our knowledge for the first time, to directly extract the overall fidelity of a heralded 
single photon by utilizing the HOM interference. 



2 Properties of the spectral correlation function 

The quantum state of photon pairs generated via PDC is usually described as 

|V>> = \0)+(J J du a dUi<Ku> a ,Wi) \u a ,ut) +0(C 2 ) , (1) 

where |£| 2 <C 1 is the pair creation probability proportional to the pump power, and the labels s,i 
refer to signal and idler. The spectral correlation function <p(oj s , Wj) arises as a consequence of energy 
and momentum conservation. The former is constrained by the spectrum of the pump pulse ot(u> s , uii) 
and the latter, being governed by the dispersion of the non-linear optical medium of the length L, is 
characterized by the phase-matching (PM) function $(w s , ui. t ) = sinc(Ak(uj s ,uj l )L/2)e iAk ^ B ' UJ ^ L ^ 2 . 
We study the spectral correlation in the vicinity of the central frequencies of pump, signal and idler w° 
(fj, = p, s, i), at which hul = Huj s + hw? and Ak° = fc s (w°) + fc 4 (w°) - fc p (w°) - x = 0. The grating 
period A of the quasi-phasematched structure can be used for adapting the desired central wavelengths. 
(See e.g. (13)). 

For simplicity, we describe the spectrum of the pump by a Gaussian envelope 

a(v s , Vl ) = e-M^)\ (2) 

where = — is the detuning from the central wavelength. In addition, we expand the phase- 
mismatch Ak(uj s , u>i) appearing in the PM function in first order Taylor series around and use the 
Gaussian approximation for the sinc-function. Thus, the PM function can be estimated by 

where k m = k' (u)®) — k' p {uJ^) is determined by the group-velocity mismatch of the non-linear medium 
and the parameter 7 = 0.193 adapts the width of the Gaussian to the width of the original sinc-function. 
Finally, the amplitude of the spectral correlation function is written as 

10(^5,^)1 = Ol(v s , Vi) \<b(v„,Vi)\ 
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which represents an ellipse in the [y s , z^)-space as illustrated in Fig. [TJa). For a decorrelated and 
separable photonic state the ellipse is oriented parallel to the (y s , ^-coordinates, i.e. the matrix M. 
in Eq. Q is diagonal. Non-vanishing off-diagonal terms correspond to a correlated photonic state. We 
stress that these properties can be tailored by modifying the pump bandwidth or the PM properties, e.g. 
by adjusting the waveguide length. 

The width of the ellipse axes are determined by the eigenvalues of Ai and are equal to 
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whereas the tilt can be extracted from the eigenvectors. In the limit of spectrally broadband pump 
(1/cr 2 0) the tilt 8 is given by tan(0) = K s / Ki and the minor axis of the ellipse is purely determined 
by the PM properties, 

1 L 2 

_ = 7 ( K 2 + K 2 ) (6) 
°PM 41 

In Fig. [TJb) we show the behavior of this width as a function of L considering the expected material 
properties of our WG. Even though this theoretical prediction gives an indication about the PM prop- 
erties, the actual spectral correlation function needs to be accurately characterized. This is especially 
important for waveguided PDC as the non-ideal properties of the guide typically influence the PM func- 
tion significantly. 
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Figure 1 : (a) Illustration of the spectral correlation following from the pump and PM properties, (b) Pre- 
dicted FWHM (full width at half maximum) of the PM at 796 nm with respect to WG length. Employed 
WG parameters: 4x6 um 2 cross section and 7.46 nm poling period. 



In terms of Fock state preparation, the joint amplitude of signal and idler in Eq. Q has to be de- 
coupled. To optimally achieve this goal we employ a short WG, thus broadening the PM bandwidth, 
and a femtosecond pump, which allows the state decorrelation via a proper design of the spectral filter. 
Filtering reduces the width of the ellipse's major axis and modifies the correlation such that the ellipse 
evolves into a circular shape or has a very small remaining tilt. Thence, the bandwidth of the filter has 
to be matched with the width of the ellipse's minor axis. 

2.1 Simple estimation of the ellipse's angle and axis 

One of the most direct ways to characterize the tilt of the spectral correlation ellipse is to measure the 
unconditional spectral marginal distributions of signal and idler. The tilt is determined by the simple 
relation tan(0) = ^j-, where Aw,, (/i = i, s) is the FWHM of the marginal. However, several spatial 
modes can propagate in the WG and an accurate determination of the tilt via this method becomes 
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difficult [33) . Another possibility relies on detuning the pump wavelength and observing the respective 
change in the central wavelengths of signal and idler. Similarly, we can also employ narrowband filters 
in signal and idler arms, and record the wavelengths at which photon-number correlation appears, thus 
mapping the spectral correlation function. 

In all PDC measurements presented in this paper we employed a 2.1 mm long, type-II, periodically 
poled KTP WG, which was pumped by a frequency doubled TkSapphire laser operating at ultrafast 
regime (796 nm, 10 nm FWHM, 170 fs autocorrelation length and 1 —4 MHz repetition rate). For 
the estimation of the ellipse's tilt we separated signal and idler photons in a polarizing beam splitter 
and launched them into single mode (SM) fibers. The single photon spectra was measured with a 
sensitive spectrograph (Andor) and the photon-number correlation with avalanche photo diodes (APDs) 
as illustrated in Fig.|2|a). First of all, we measured the spectral marginal distributions of signal and idler 
as shown in Fig. |2|b). After that, two narrowband interference filters were set on the beam paths of 
signal and idler and the photon-number correlation was mapped by tuning the filters. For the accurate 
determination of the ellipse's tilt the filter wavelengths were recorded with the spectrograph. From the 
results shown in Fig. |2jc) we extracted the following value for the tilt: 9 = 54.7° ± 1.5°. 
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Figure 2: (a) Setup for measuring the tilt of the PM function. [The list of abbreviations appears in Fig. [5]] 
(b) Nearly degenerate, measured marginal distributions of signal and idler and (c) the wavelengths at 
which coincidence counts between signal and idler were recorded. 



Regarding the major axis of the correlation ellipse, note that its width is highly dependent on the 
bandwidth of the pump as illustrated in Fig.[TJa). The FWHM of the major axis can be roughly estimated 
from the marginal distributions by Aui s / cos(8) w Aw, / sin(0). Our results indicate that pumping 
the WG with 1 — 1.5 nm broad pulse leads to an approximately 35 — 40 nm broad major axis, which 
increases to values between 55 — 60 nm, if the pump is 2 — 2.5 nm broad. 

According to Fig.[TJb) the exact determination of the PM bandwidth of our WG, i.e. the minor axis 
of the ellipse, requires sub-nanometer resolution. Instead of employing sensitive monochromators, we 
studied the second harmonic (SH) response of the WG. In this process two photons at the fundamental 
frequency give rise to a single up-converted one. Therefore, the fundamental probe always lies on the 
line Wj = uj s of the (lo s , (jj)-space, as illustrated in Fig.[5Ja). In other words, tuning the frequency of 
the fundamental probe and measuring the respective SH response allows us to characterize the width of 
the PM function. If the fundamental probe is narrow with respect to the PM bandwidth apui an d the 
tilt of the PM only slightly deviates from 45°, the envelope of the SH response is approximately equal 
to 2 

I(i*h) ~ exp f-^V (7) 

We recorded the SH response with a grating spectrometer by tuning the fundamental probe (795.7 nm, 
0.6 nm FWHM, and 2.6 ps autocorrelation length) in steps of approximately 0.2 nm over the phase 
matching curve [setup shown in Fig.|3jb)]. In our case, the measurement is limited by both the res- 
olution of the spectrograph (0.36 nm at 800 nm) and the bandwidth of the fundamental probe beam. 
Nevertheless, our results shown in Fig.|3jc) indicate that we can probe the shape of the PM function and 
that its bandwidth is on the order of the probe bandwidth. 
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Figure 3: (a) Illustration of how to scan the PM structure with SH measurements, (b) Setup employed 
in the SH measurement. [The list of abbreviations appears in Fig. [5]] (c) Measured SH responses to 
different fundamental wavelengths. 



By combining our findings we are able to reconstruct the correlation ellipse \4>{v s , Vi)\ 2 of our pho- 
ton source. Using the estimated values of the tilt and width of the PM, and taking into account a pump 
with FWHM of 2.5 nm we obtain the ellipse shown in Fig.Qa). The result implies that signal and idler 
are highly frequency correlated, which is a consequence of the long interaction length of the PDC pho- 
tons in the WG. Further, our result suggests that the optimal filter bandwidth is on the order of 0.6 nm. 
Unfortunately, spectral filters at 800 nm with high transmission and narrow bandwidth are rarely avail- 
able. In consequence, we study in the next section the decorrelation of the generated photonic state by 
employing spectral filters with 1 nm and 2.5 nm bandwidths. We expect slight remaining correlations 
after sending the prepared state through these filters as shown in FigsQb) and|4|c), respectively. 
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Figure 4: (a) The reconstructed spectral correlation ellipse \4>{v s , Vi)\ 2 - We applied in the simulation 
2.5 nm broad pump at centered at 795.7 nm. Correlations remaining between signal and idler after (b) 
1 nm and (c) 2.5 nm broad spectral filtering. 



3 Decorrelation of signal and idler 

The indistinguishability between signal and idler can be studied via a HOM interference, where signal 
and idler are superimposed in a 50/50 beam splitter. The full visibility is achieved either if the correlation 
function is circular, or if it is symmetric with respect to the axis of energy preservation, i.e. if the 
correlation ellipse is oriented at 45° [34) . Here we employ this measurement to test and corroborate the 
acquired knowledge about the spectral ellipse. 

We adapt the analysis of Avenhaus et al. 1 12 1 to the study of the HOM visibility with spectral filter. 
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The interference between signal and idler can be characterized by the overlap factor 



O 



J J dv s di/i<l>(u s ,Vi)(l>*(vi,u s ) 
J J dv s dvi \§(y s ,Vi)\ 



(8) 



Using the Gaussian approximation of the spectral correlation function, the overlap factor can be ex- 
pressed as a product of a spectral (s) and temporal (i) overlap i.e., O = O s Ot- While spectral mismatch 
results from an asymmetric spectral correlation function, temporal mismatch originates from the phase 
terms, which are ultimately related to the properties of the birefringent media. These can be inferred 
from the properties of the correlation ellipse via Eq. pj. If perfect overlap was necessary, the temporal 
mismatch could be compensated with an optical delay line. 

The defined overlap factors depend on the ellipse's tilt 6 and aspect ratio A, defined as the ratio 
between the major and minor axes. They can be expressed as 
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yj(l + .4 4 )(1 - sin 2 (26)) + 2A 2 (l + sin 2 (26)) 



and 



(9) 



O t = exp 
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A 2 [l-sin 2 (2fl)]„4 2 + [l-sin(2fl)] 2 
'2q [1 + A 4 } [1 - sin 2 (26)} + 2 A 2 [1 + sin 2 (20)] , 
where we take use of an estimation that the PM width and the ellipse's minor axis are equal. 

Since the twin-photons of a type-II process have orthogonal polarizations, we measured HOM in- 
terference fringes between signal and idler by rotating a half-wave plate in front of a polarizing beam 
splitter and recording the coincidence counts between the two output arms. The experiment was realized 
in three different cases: (I) without spectral filter, (II) including in the setup a spectral filter with FWHM 
of 2.5 nm and, (III) decreasing the width of the filter down to 1 nm. The filters were inserted in the joint 
beam path of signal and idler as shown in Fig. |5ja). According to auxiliary experimental observations, 
the time delay between p- and s-polarized modes is on the order of 700 fs. Since the temporal width 
of signal and idler, estimated by the filter bandwidth, is on the order of 1 — 2 ps and no delay line is 
inserted in the setup, the temporal overlap cannot be neglected in our measurements. 
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Figure 5: Setup for preparing and characterizing single photons from waveguided PDC source applied 
for (a) measuring HOM interference between signal and idler, (b) recording the statistics of the heralded 
state, and (c) measuring the HOM interference between two independent sources. Abbreviations: SHG 
= second harmonic generation (0.5 mm thick BBO), L = lens, P = pinhole, DM = dicroich mirror, CF 
= color glass filter, HWP = half-wave plate, PBS = polarizing beam splitter, WG = waveguide, IF = 
interference filter, ND = neutral density filter, FC-PBS = fiber coupled PBS, APD = avalanche photo 
diode, SG = spectrograph* and SM = spectrometer*. (* only shown in Figs|2]and[3]) 



The visibility is defined by V = (C max — C m i n ) / (C ma x + C m i n ), where C indicates the amount 
of coincident counts. In terms of the overlap factor, the visibility can also be expressed as V — 
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(1 + 0)/(3 — O). The classical case, in which signal and idler are completely distinguishable, cor- 
responds to the situation of zero overlap, i.e. V = 1/3. We measured, without any background subtrac- 
tion, the following visibilities: V(I) = 0.34 ± 0.02, V(II) = 0.58 ± 0.02, and V(III) = 0.81 ± 0.03, 
corresponding to each of the above mentioned cases. The corresponding overlap factors must be com- 
pared with the theoretical prediction given in Eqs d9l) and ( 10 1. For this we approximate the aspect ratios 



in a straightforward manner. When no filtering is applied, A is estimated as the ratio of the experimen- 
tally extracted major and minor axes, i.e, A(I) ~ 90 — 100. When spectral filtering is applied, we 
determine A as the ratio of the filter bandwidth to the PM width i.e., A{II) — 4.2 and A(III) — 1.7. 
As expected, the spectral overlap alone does not explain our results, as can be seen in Fig. [6] By includ- 
ing the temporal overlap in the analysis, the employed model predicts the right tendency, even though 
in the experiment we do not have any compensation for dispersion or chirp. 
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Figure 6: Measured overlap values (squares) together with theoretical estimation for the spectral (dashed 
line) and total (solid line) overlaps with respect to the aspect ratio. We consider the PM tilt of 9 = 
54.7° ± 1.5° (indicated with the shaded areas). 



We conclude that our model, which includes only experimentally extracted characteristics of the 
PM, explains with high confidence results of an independent experiment. In addition, it predicts how 
the spectral correlation between signal and idler can be manipulated with the selected filter. In the next 
sections we study how the heralded state is affected by the higher photon-number contributions and 
remaining spectral correlation. 



4 Measurement of the heralded statistics 

In the low gain regime heralding of pure signal states is normally realized by considering a single 
spectral filter at the trigger. However, at higher gains multiple modes are excited simultaneously. Thus, 
filtering only the trigger results in the detection of uncorrelated signal modes. In addition to this, the 
higher photon number contributions decrease the fidelity of the heralded single-photon state. We studied 
these effects by measuring, how the heralded statistics depends on whether or not a narrowband spectral 
filter is added in the signal arm. As before, we pumped our PDC source with a frequency doubled 
TkSapphire laser. The pump beam was sent through a simple spatial mode cleaner before being launched 
into the WG. After that, the pump was filtered out and the twins were divided into two independent paths. 
The idler was sent through 1 nm broad interference filter and, finally, coupled into SM fiber. 

We measured the heralded statistics in two different configurations: either without filter or with a 
2.5 nm interference filter inserted in the signal arm. The setup shown in Fig. |5jb) was arranged such 
that the idler was used as a trigger and the statistics of the signal was measured with a "two-bin TMD". 
This detector consists of a single 50/50 beam splitter, constructed by using a half-wave plate in front of 
a polarizing beam splitter. The loss tolerance of the TMD was achieved by accurate calibration of loss 
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Figure 7: Heralded statistics at different pump powers with (a) no filter and (b) 2.5 nm filter in signal, 
(c) Mean photon number of the statistics with respect to the pump power: diamonds referring to plot (a) 
and squares to plot (b) are shown together with linear fits to the data. 



L{rf) and convolution matrices C [35]. The signal statistics g was then deduced by inverting the click 
statistics pdick = CL(rj)g with a maximum-likekihood technique ]36) . 

From our results one can observe that the heralded statistics presents a large two-photon contribution 
[Fig.[7|a)], which is reduced by placing a narrow band filter in the signal arm [Fig.[7Jb)]. We understand 
this result as an indication that, at higher gains, a large amount of distinct spectral modes are excited 
simultaneously. If two idler photons are generated there is a high probability that just one of them has a 
good overlap with the filter mode. Therefore, by using a filter in the signal arm it is possible to reject the 
signal photon created together with the spectrally mismatched idler. The click statistics were inverted 
using the calibrated efficiencies, which in both cases were 9%, corresponding to the ratio of coincident 
counts to singles. The highest trigger rate was 5.9 kHz, which is equivalent to a heralding probability of 
0.16%. Regarding the single-photon state preparation, the one-photon contribution of the filtered signal 
was 95% at this trigger rate. 

A comparison between the filtered and unfiltered statistics allows us to extract the degree of mode 
reduction given by the filter. For that, we model the multimode statistics as the convolution of several 
uniformly occupied spectral modes [16,37]. We begin the analysis with the single-mode case. The 
heralded statistics, being conditioned on a trigger event, has the form 

Pi PnHclick / 1 1 \ 

n\click = —j^ ■ (ll) 

^click 

Moreover, the joint probability of having n photons in signal and detecting a click from idler is given 
by Pnnciick = P„ [1 — (1 — where rjt is the efficiency at which the filtered mode is detected 

and P„ the probability of creating n signal (idler) photons, which always occur in pairs for PDC. 
Since in the experimental realization losses are very high, the joint probability can be approximated by 
Pnndick ~ r] t nP£ DC . The overall click probability is estimated as P cHck = Y^ n P nnclick ~ Vt (n), 
where (n) indicates the mean photon number of the marginal PDC statistics. In this loss regime, the 
heralded statistics become independent of the trigger efficiency rjt and its mean photon number can be 
evaluated as (riheralded) — ^2 n n Pn\ciick — (n 2 ) I ( n ), with (n 2 ) being the second moment of the 
marginal PDC statistics. 

In the case of a multimode PDC with several uniformly excited modes, the statistics P^ D is given 
by the convolution of the modes, each having the mean photon number (n sm ) and the second moment 
(n 2 sm ). For M excited modes we find that (n) = M (n sm ) and (n 2 ) = M (n 2 am )+M(M -1) (n sm ) 2 , 
from which we can calculate the multimode mean photon number of the heralded statistics. Considering 
a thermal single mode marginal distribution, g^ 71 = |x| 2 "(l — \x\ 2 ) with |x| being the gain of the PDC 
mode, the mean photon number of the heralded statistics can be approximated by 

/ n 2 \ 

("heralded) = / ~\ + ( M - 1) (>W) 



8 



« 1 + (M +l)|x| 2 ■ (12) 

We conclude that, in the limit of low pump powers, such that |x| 2 oc pump power, the slope of the mean 
photon number depends on the number of modes. In Fig.|7|c) we show the mean photon number of the 
measured heralded statistics with respect to the pump power. By comparing the slopes, we estimate that 
the reduction in the number of modes due to the filter at signal is on the order of 30. We emphasize that 
the model used is an approximation, where all modes are treated similarly, even though the trigger mode 
is a special one. 



5 Interference between independent sources 

In this section we study the HOM interference between two independent sources: the heralded signal 
and a strongly attenuated coherent field, here called reference. The visibility of the interference dip 
depends, on the one hand, on the presence of higher photon-number components and, on the other 
hand, on the spectral overlap between the target state and the reference. We experimentally investigate 
the cases of two- and three-fold HOM interferences. The former case corresponds to the interference 
between signal and reference neglecting the PDC photon-number correlation, while in the latter case the 
signal state is heralded by trigger events. 

We begin our study by considering a spectrally pure signal state. Further, we decompose the signal 
only up to the second order in the photon-number basis, which is justified under the experimental con- 
ditions. After that, we complete our description by adding the effect of spectral correlations between 
signal and idler. In the experiment, signal and reference are launched into the ports a and b of a 50/50 
beam splitter. The density matrix of the input state p a b can be written as 

Pab = (po |0> <0| + Pl |1> (1| + P 2 |2> <2| ) a ® ( |/3) (f3\ ) & , (13) 

in which p n (n = 0,1,2) corresponds to the statistics of the photon-number mixed signal mode 
with J2 n P n = 1- The re f erence field has amplitude j3 and spectral distribution u(ui), such that 
fdu> \f3u(uj)\ 2 — \{3\ 2 - Since mode "a" is spectrally pure, its one- and two-photon components can be re- 
spectively described by | l) a = fdujf(uj) a\uj) |0) o and|2) o = ^ Jdu) J du> f(w) f(£j) a\uj) d\Cb) |0) a , 
in which the spectral distribution f(uj) is normalized, i.e. J dcu \f(oj)\ 2 = 1. 

Interference does not take place when the two input pulses arrive at the beam splitter at different 
times. To model that, we explicitly write the evolution of the fields considering a time delay r. In the 
frequency domain, the beam splitter relations are given by 

at (a;) = -L[ctH-dt M 

St (a;) = i= [et( w ) +( ?t( w )] , (14) 

in which c and d are the two output modes. These equations allow us to calculate the evolved state 
p c d- Finally, the probability of a coincident click can be calculated via P = Tr |/5 C( i fl c <g) ft<j|, where 

f[ c ® fl d = (1 - |0) CC (Q|) ® (1 - |0)d rf(0|) describes the POVM of two simultaneous clicks in the 
APDs. Using the bosonic commutation relation [a(ui), a'(ui)] = 5(ui — tl>) and the fact that b(oj) \ f3) = 
(3u(u>) |/3), the probability of detecting a coincident click is given by 

P = Po(l-e-^ 2 ) 2 +Ml-e-^ 2 )- Pl T(r)\f3'\ 2 e-^ 2 

+ P2 (i - 2 eH/5 '' 2 ) - P 2 eH/5 '' 2 ( T ^'\ 2 + T '( T ) } ' (15) 

with /3' = P/y/2. The functions T(t) and T'(r) have the form 

T(t) = Jdujx Jdu 2 u*{uj x ) g(ui,cj 2 ) u{io 2 )e lT{ ^- U2) , (16) 

T'(t) = (dux / du> 2 I duj 3 I duj 4 u* (u} X )u* (w 2 ) K^, u 2 , u>a, u>4.) u{uj 3 ,)u{Lu i )e lT{u>1+UJ2 - UJ3 - UJi \ 
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Figure 8: Contributions affecting the HOM interference between signal and reference. 



where g(uo u oj 2 ) = f(ui)f*(u) 2 ) and h(ui,u 2 ,u 3 , W4) = f (^1) f (^2) f* (^3) f* (^4) ■ At origin (r = 
0) these factors determine the spectral overlap of the one- and two-photon wave-packets with respect to 
the reference. Finally, the visibility of the interference dip is defined by 



V 



P(t -> 00) - P(t = 0) 
P(t 00) 



(17) 



In order to simplify our model we compare the independent impact of the one- and two-photon com- 
ponents of signal on the visibility. Considering perfect spectral overlap between these components and 
the reference field, and appliying values of |/3| 2 compatible with the experimental realization, we esti- 
mate that the interference of the one-photon component would have a visibility of above 98%, whereas 
the one of the two-photon component would be less than 5%. Therefore, in the regime of power consid- 
ered here, the two-photon component of signal gives only a background contribution to the coincidences. 
Thus, we consider the spectral degree of freedom exclusively for the one-photon component of signal 
and reference. Further, we approximate e~l^ I » 1 — |/3'| 2 , which is well valid within the experimental 
frame, and neglect the terms arriving from the interference of more than two photons (total number in 
mode "a" plus "6"). The remaining contributions are shown in Fig. [H] and they can be retrieved from 
Eq. ( fT3] > as follows. The first term corresponds to accidental counts due to the two-photon component of 
the reference and it can be re-expressed as po(l — e - '' 3 ' ) 2 po|/3'| 4 [Fig.^a)]. The second and third 
terms account for the events where both signal and reference contain one single photon [Fig.[HJb)]. The 
probability of simultaneous accidental clicks due to one photon from signal and another from reference 



is given in the second term as p%(l — 



-|/3'| : 



) ~ Pi I/?' I ■ Considering the first-order photon-number 



contribution, the third term describing the quantum interference between the one-photon components 
of signal and reference can be rewritten as piT(t)\(3'\ 2 . The fourth term takes into account acciden- 
tal counts coming from signal, p 2 (l ~ 



-l/3'l 



) w ip 2 [Fig. [sjc)] . The last term, corresponding to 



interference of at least three photons, is an order of magnitude smaller the other terms and thence it is 
neglected. 

With the described simplifications, Eq. ( 15 1 can be rewritten as 



1 



P(r) =Po|/3? - T{r)) + -p 2 



(18) 



After substituting Eq. ( p"8j ) into Eq. ( 17 1, the visibility can be written with respect to the mean photon 



number of the reference 



(19) 



which reaches its maximum at 1/31^ = ^/2 p 2 /pa- 

In the experimental situation the joint-spectral distribution of signal and idler is not separable, and 
thus the prepared state is spectrally mixed. Our model can be adapted to this condition by replacing the 
function g(uj\ , W2) by a new one, in which a trace operation over the idler mode is realized. The spectral 
density of the PDC state after filtering is given by 



1 

N 



(20) 
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Po(%) 


Pi(%) 


P2(%) 


2-fold 


99.7896 


0.2101 


0.0003 


3-fold 


94.920 


5.065 


0.015 



Table 1: Measured signal state statistics p n (n = 0, 1, 2) without loss inversion. 



eo(%) 




ea(%) 


3-fold 0.2±0.2 


93.1±0.6 


6.6±0.6 



Table 2: Statistics g„ of the heralded state in Table [T] inverted with 4.8% ± 0.2% efficiency. 



where the t^ioS) is the intensity transmission profile of the signal and idler filters employed in the state 
generation and N accounts for the normalization. 

Apart from the visibility, the width of the HOM interference dip offers information about the spectral 
properties of the prepared state. Using the Gaussian approximation for the spectral correlation function, 



reference, and filter profiles, the width of the HOM interference dip can be estimated via Eq. (18i 
Disregarding the background contributions, it can be written as 



Spectral (r) (X 1 - T(t) = 1 - T max expl-^1 , (21) 
where the temporal width of the interference dip, ai = \ + \ H — \ — h sl " 6 , is determined by the 

° ° '/3 a t 3 a pm 

inverse widths of the pump a, reference ap, signal filter cr ts , and phase matching a pm , as well as by the 
ellipse's tilt 9. The amplitude T max describes the maximal expected value of the overlap, ranging from 
zero to one. We emphasize that the two- and three-fold interference dips have the same temporal width. 

For the measurement of the HOM interference we used the setup shown in Fig. |5jc). A portion 
of the original laser beam was heavily attenuated and launched through a polarization maintaining SM 
fiber-coupler in s-polarization. The signal was launched to the fiber in p-polarization and overlapped 
with the reference at a 50/50 beam splitter, consisting of a half-wave plate and a polarizing beam splitter. 
Two 1 nm broad interference filters were inserted in the setup: one in the path of the trigger and the other 
in the joint arm of the signal and reference. We measured the HOM interference dip in the two- and 
three-fold cases for different amplitudes of the reference. By blocking the reference arm we measured 
the statistics of the signal state and obtained the results summarized in Table [T] 



Substituting the measured statistics back into Eq. ( 19 1 we can find the maximal expected visibilities 
of V™J ld = 0.46 and V™*^ = 0.75. In these values we disregard the spectral mismatch, i.e. T = 1, 
and consider only the degradation due to the two-photon component. Our experimental results are 
shown in Fig. [9] where we plot the measured visibilities with respect to the mean photon number of 
the reference beam \(3\ 2 . The latter can be found by measuring the single-counts rate of the reference 
in one of the detectors while blocking the signal, P S i ng ios = 1 — <s !/ 2 1 ^1 « |/3| 2 /2. The discrepancy 
between the expected and observed visibilities is a consequence of the spectral mismatch. We extract 



the value of the spectral overlap T from the measurement by fitting Eq. ( 19 1 to the experimental data. 
The obtained overlap values are T2-fdd = 0.41 ± 0.02 and T3_f id = 0.65 ± 0.03. This clearly shows 
the increase in purity achieved via filtered herald, i.e the degree of spectral correlation between the 
twins is reduced for the selected conditioned events. Finally, the fidelity of the prepared state is given 
by T = J (l|p|l) = y/Tgi, where T is the spectral overlap and Qi is the one-photon component of 
the inverted statistics shown in Table [2] Combining the results of both measurements we find that the 
fidelity of the heralded state was J-3-foid = 0.78 ± 0.03. 

In FigsfTOfa) andfTOfo) we show example measurements of the two- and three-fold interference dips 
that were recorded with reference powers close to the optimized values |/3|^ ax in Fig. [9] The measured 
visibilities were Va-foid = 0.18 ± 0.02 and VWoid = 0.48 ± 0.03, respectively. The temporal width of the 
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0°00 0.01 0.02 0.03 0.04 0.05 
Mean photon number of reference, |/?j 2 



Figure 9: Visibility of the HOM interference dip with respect to the mean photon number of reference. 
Measured points are marked by diamonds (two-fold interference) and squares (three-fold interference). 
Solid and dashed lines correspond to the fitted curves from which the spectral overlap value can be 
extracted. 



interference dips was 2.0 ± 0.2 ps. Following Eq. pT| ) we theoretically study how this width changes 
as a function of the PM bandwidth [Fig. [TOfc)]. The comparison of the measurement with our model 
suggests a FWHM of 0.5 ± 0.1 nm for the PM, which is in good agreement with the value extracted 
from SH measurements. Further exploring our knowledge about the spectral correlation function, we 
use Eqs. ( 16 1 and ( |20| in order to estimate the maximal spectral overlap T max in the two- and three- 
fold cases. In Fig.|10[c) we illustrate the theoretical predictions and compare them with the measured 
overlap values. In the three-fold case we find a good agreement with the model. However, the measured 
value of the two-fold overlap deviates largely from the expected one. As no selection of the detection 
events by triggering is applied in the two-fold case, we believe that the properties of non-ideal WG, such 
as fluorescence and higher order spatial modes with complex spectral structure, degrade the result. In 
addition to this, the expected values have not been achieved due to other experimental imperfections, in 
particular mechanical instabilities and spatial mode mismatch. 
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Figure 10: Measured HOM interference dips between the prepared signal and the reference in (a) two- 
and (b) three-fold cases, (c) Theoretical estimation of the temporal width (black solid line) and maximal 
spectral overlaps (three-fold: dashed line, two-fold: dotted line) with respect to the PM bandwidth 
(employed FWHM bandwidths: pump 2.5 nm, reference and filters 1 nm, and PM tilt: 9 = 54.7°). 
Horizontal lines refer to the measured values of overlap. Shaded area indicates the FWHM of the PM 
estimated by the temporal dip widths. 
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6 Conclusions 



We have studied the ability to tailor the spectral properties of a conventional PP-KTP waveguide. Infor- 
mation about the spectral correlation between signal and idler was retrieved from simple experiments, 
enabling to model our system with a standard theory of parametric down-conversion. The viability of 
this model was tested via the Hong-Ou-Mandel interference between the twin beams. Since the produc- 
tion of more evolved non-Gaussian states, as displaced single-photon Fock states, requires an overlap 
between the signal and an independent reference field, we studied the quantum interference between the 
heralded state and a coherent reference field. The spectral overlap was directly extracted from the vis- 
ibility of the Hong-Ou-Mandel interference. Further, we experimentally characterized the statistics of 
the heralded state. Considering both the spectral overlap and the one-photon contribution, we obtained 
by a direct detection the fidelity of 78% between the prepared state and a single-photon Fock state. In 
addition, our observations show that heralding decreases the number of modes only if both signal and 
idler are spectrally filtered. Although we found that in the current experiment the spectral overlap be- 
tween the reference and the signal state is not yet optimized, we provide a clear and simple recipe for 
tailoring the characteristics of the spectral correlation and for testing the degree of decoupling, i.e., the 
fidelity of the heralded state. 
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